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ABSTRACT 

The central regions of the gravitationally lensed quasar Q2237+0305 can be indi- 
rectly resolved on nano-arcsecond scales if viewed spectrophotometrically during a 
microlensing high magnification event (HME). Q2237+0305 is currently being moni- 
tored from the ground (eg. OGLE collaboration, Apache Point Observatory), with the 
goal, among others, of triggering ground and spacecraft based target of opportunity 
(TOO) observations of such an HME. In this work we investigate the rate of change 
(trigger) in image brightness that signals an imminent HME and importantly, the 
separation between the trigger and the event peak. In addition, we produce colour 
dependent model light-curves by combining high-resolution microlensing simulations 
with a realistic model for a thermal accretion disc source. We make hypothetical target 
of opportunity spectroscopic observations using our determination of the appropriate 
trigger as a guide. We find that if the source spectrum varies with source radius, a 
3 observation TOO program should be able to observe a microlensing change in the 
continuum slope following a light-curve trigger with a success rate of > 80%. 

Key words: gravitational lensing - microlensing - numerical methods - accretion-disc. 



1 INTRODUCTION 

The QSO 2237-1-0305, sometimes known as Huchra's lens 
(Huchra at al. 1985), is perhaps the most remarkable gravita- 
tional lens yet discovered. It comprises a foreground barred 
Sb galaxy (z=0.0394) whose nucleus is surrounded by four 
images of a radio-quiet QSO (z=1.695). Ground based spec- 
troscopic observations have verified that all four images are 
similar QSO's at the same redshift (Adam et al. 1989). 
These images are separated by 0.7-0.9" from each other and 
the nucleus. Broad band monitoring has shown that signif- 
icant microlensing events occur (eg. Irwin et al. 1989; Cor- 
rigan et al. 1991). Since the optical depth to microlensing 
is of order unity at each of the image positions (eg. Kent 
& Faico 1988; Schneider et al. 1988; Schmidt, Webster & 
Lewis 1998), the magnification effects on the source can be 
considered as a network of caustics moving across the source 
plane. Q2237-I-0305 provides a unique opportunity to use mi- 
crolensing events to study the geometrical structure of the 
background QSO for two good reasons. Firstly, the relative 
closeness of the lensing galaxy means that the time delay 
between the four quasar images is less than a day. Thus it is 
easy to separate microlensing events from intrinsic variations 
of the QSO. Also, for this QSO the time-scale of microlens- 



ing events is typically 30-50 days, reduced by a factor of 10 
from typical time-scales of perhaps a year. 

Detailed modelling (eg. Wambsganss & Paczynski 1991; 
Fluke & Webster 1999) shows that microlensing events re- 
sulting from the source crossing a caustic will show spectral 
variations which depend on the relative sizes of different 
emission regions at different wavelengths. Ground based ob- 
servations have confirmed differential amplification of the 
emission region. Lewis et al. (1998) determine the ratios of 
emission line equivalent widths relative to one image. They 
show (i) the ratios remain fairly constant for one image from 
line to line, suggesting that the sizes of the emission regions 
for the lines are not greatly different, (ii) that the ratios 
vary from image to image for a single epoch by a factor of 
~ 2.5, and (Hi) that the ratio for a single image varies as 
a function of time, i.e. as a result of a microlensing event. 
These results are consistent with earlier results of Fillipenko 
(1989) who measured a ~ 25% difference in the width of the 
Mgll lines between the A and B images. 

A QSO is generally described as a continuum source sur- 
rounded by small gas clouds which are photo-ionised by the 
central source. The size of the central continuum source in 
a QSO is expected to be at least ten times the gravitational 
radius (r^) of the black hole thought to power the source: 
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Tg ~ 10^^(Mi3///10** Mq) cm. This size is corroborated by 
the length scale, ct ~ IQ^'^thours cm, associated with X-ray 
and optical changes in the continuum emission. 

Reverberation mapping provides one technique for de- 
termining the geometric structure of a QSO, although the 
technique has mostly been applied to nearby Seyferts. Sub- 
stantial monitoring programs have already confirmed a pic- 
ture where the continuum region and the region emitting 
the highest ionisation lines vary faster and are therefore 
smaller than the resolution scale of the reverberation tech- 
nique. Typical sampling time-scales are several days, which 
corresponds to 2.5 x 10^^ tdays cm. 

Microlensing allows us to probe length scales which are 
an order of magnitude smaller than this, and which will not 
be resolvable by any other technique in the foreseeable fu- 
ture. Already variability amplitude (Wambsganss, Paczyn- 
ski & Schneider 1990) and time (Wyithe, Webster, Turner 
& Mortlock 2000 (hereafter WWTMOO)) constraints from 
the observed microlensing event limit the R-band contin- 

uum emission region to be less than 2 x 10^^ hjJ' cm. This 
size limit has been used by Ranch and Blandford (1991) 
to argue that the continuum emission cannot result from 
a thermal accretion disc, although Jarozynski et al. (1992) 
dispute this result. Spectrophotometric observations dur- 
ing a HME could resolve the issue. Although simulation of 
caustic crossings using realistic accretion disk models is the 
most straightforward approach theoretically and the least 
demanding of observational data quality, it is an intrinsi- 
cally model dependent mode of investigation. A difi'erent 
approach is to invert the light curve observed to directly ob- 
tain a surface brightness profile of the source (eg. Grieger et 
al. 1998; Grieger et al. 1991; Agol & Krolik 1999; Mineshige 
& Yonehara 1999). 

We discuss microlensing theory related to a proposed 
experiment to measure spectroscopic change during an 
HME. The experiment is a target of opportunity program 
with Hubble Space Telescope, where the observations would 
be triggered by the observation of the onset of a caustic 
crossing from ground based monitoring of the QSO. The goal 
is to observe a caustic passing across the continuum region, 
magnifying difi^erent length-scales as a function of time. The 
wave-lengths chosen will cover the Ly Limit, Ly/3/OVI, Lya, 
NV 01, SilV/OIV, CIV, HeH and MgH emission lines and 
a number of associated absorption lines, covering more than 
3700^ in the rest frame of the QSO. The long wavelength 
baseline will enable determination of whether the contin- 
uum emission at all wavelengths is emitted from the same 
geometric regions or whether some regions are bluer. 

This paper discusses the light-curve derivatives (trig- 
gers) in monitoring data that signal an imminent HME. 
The aim is to determine a reliable function to predicting 
the onset of an HME from monitoring data in real time. 
We then discuss one potential outcome of the TOO exper- 
iment; the spectral change of a thermal continuum during 
an HME. Sample triggering scenarios are applied to model 
microlensed light curves produced for a thermal accretion 
disc, and the probability of success of an example target of 
opportunity program is discussed. 



2 THE MICROLENSING MODEL 

Throughout the paper, standard notation for gravitational 
lensing is used. The Einstein radius of a IMq star in the 
source plane is denoted by ?7o- The normalised shear is de- 
noted by 7, and the convergence or optical depth by k. The 
basic model for microlensing at high optical depth was first 
described by Kayser, Refsdal & Stabell (1986). It comprised 
a disc of point masses having a size such that a large frac- 
tion (> 99%) of macroimage flux is recovered (Katz, Bal- 
bus & Paczynski 1986; Lewis & Irwin 1995). To construct 
a microlensed light-curve we use the contouring technique 
of Lewis et al. (1993) and Witt (1993). For the microlens- 
ing models of Q2237-(-0305 presented in the current work we 
assume the macro-parameters for the lensing galaxy calcu- 
lated by Schmidt, Webster & Lewis (1998). The standard 
notation introduced by Yee (1988) is used to describe these 
images. Where required a cosmology having 57 = 1 with 
Ho ~ 75 km sec~^ is assumed. 

We describe the microlensing rate using the effective 
transverse velocity which is defined as the transverse ve- 
locity that produces a microlensing rate from a static model 
equal to that of the observed light-curve (Wyithe, Webster & 
Turner 1999 (hereafter WWT99)). The effective transverse 
velocity therefore describes the microlensing rate due to the 
combination of the effects of a galactic transverse velocity 
and proper motion of microlenses. Our calculation of the 
light-curve derivatives that precede an HME assume that the 
effective transverse velocity accurately describes not only 
the distribution of light-curve derivatives during an HME 
(Wyithe, Webster & Turner 2000a (hereafter WWTOOa), 
but also the distribution of orientations between the caustic 
and source trajectory, and hence the event duration. 

We have previously obtained the following normalised 
probability distributions. These were obtained under the as- 
sumption that the source size S <^ rjo. Evidence in favour of 
this assumption was presented in Wyithe, Webster & Turner 
(2000c). 

i) Pa(S'l(m), Ue//), the probability that the continuum source 
diameter is between S and S + dS given a mean microlens 
mass (m) and an effective galactic transverse velocity v^ff 
(WWTMOO). 

ii) Pv{veff\{m)) the probability that the effective galactic 
transverse velocity is between Vef / and Vef / -I- dvef / given a 
mean microlens mass (m) (WWT99). 

Hi) Pm{{m)), the probability that the mean microlens mass 
is between (m) and (m) +d{m) (Wyithe, Webster & Turner 
2000b). 

Pvi^ef f\{m}) and pm{{m)) were computed using flat 

{p{Vtran) <X dVtran), and logarithmic {p{Vtran) oc " ) 

assumptions for the Bayesian prior for galactic transverse 
velocity (Vtran). Pv{vef f\{m)) was found to be insensitive 
to the prior assumed, however Pm{{m)) showed some de- 
pendence. In the remainder of this paper we use pm{{m)) 
calculated using the assumption of a logarithmic prior. We 
note that the assumption of the flat prior raises the average 
light-curve derivative by a few percent. 

The functions ps{S\{m),Veff), Pv{veff\{m)), pm{{m)) 
and the HME statistics presented in this paper were com- 
puted for the following assumptions of smooth matter den- 
sity, photometric error, and direction of the galactic trans- 
verse velocity. Two models are considered for the distribu- 
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tion of microlenses, one with no continuously distributed 
matter, and one where smooth matter contributes 50% of 
the surface mass density. Two orientations were chosen for 
the transverse velocity with respect to the galaxy, having 
source trajectories parallel to the A— B and C— D axes. The 
two orientations bracket the range of possibilities, and be- 
cause the images are positioned approximately orthogonally 
with respect to the galactic centre correspond to shear values 
of 7a,7b < 0,7c, 70 > and 7A,7fl < 0, 70,70 > re- 
spectively. Photometric error was simulated by perturbing 
the model light-curve by an amount distributed randomly 
from a Gaussian of halfwidth a. The simulations used two 
different estimates of the error in the photometric magni- 
tudes. In the first case a small error was assumed (SE). For 
images A and B, (7se=0.01 mag, and for images C and D 
o"S£;=0.02 mag. In the second case, a larger error was as- 
sumed (LE). For images A and B, (t_l_e=0.02 mag and for 
images C and D o"_l_b=0.04 mag. The random component of 
observational error quoted in Irwin et al. (1989) was 0.02 
mag. 

Both the microlensing rate due to a transverse velocity 
(eg. Witt, Kayser & Refsdal 1993), as well as the corre- 
sponding rate due to proper motions (WWTOOa) are not 
functions of the details of the microlens mass distribution, 
but rather are only dependent on the mean microlens mass. 
We therefore limit our attention to models in which all the 
microlenses have the same mass since the results obtained 
will be applicable to other models with different forms for 
the mass function. 



3 INTRINSIC VARIATION OF Q2237+0305 

Intrinsic variation of the quasar in the 22374-0305 system 
cannot be directly measured since one cannot be sure how 
much variation is due to microlensing and how much is 
intrinsic to the source. The afore-mentioned analyses of 
microlensing in Q2237-I-0305 to determine ps, Pv and Pm 
have utilised difference light-curves. This is practical for 
Q2237+0305 since the relative time-delays between images 
(< 1-day (eg Schneider et al. 1988)) are small compared with 
typical sampling time-scales (weeks or months). However our 
current goal is to determine whether a particular observed 
light-curve derivative signals a forthcoming HME. Since an 
observed (single image) light-curve derivative depends on 
the super-position of microlensing and intrinsic source varia- 
tion, the calculation of an event trigger requires the inclusion 
of a statistical description of the intrinsic source variation. 

Giveon et al. (1999) analysed optical variability in the 
sample of Palomar-Green quasars. They find that over time- 
scales between 100 and 1000 days the sources have a power 
spectra that is of a power-law form. Pi, oc f'" with 71 ~ 
—2. We model intrinsic variability for Q2237-I-0305 using a 
power-spectrum of power-law form. Since long duration in- 
trinsic variability may have been observed in Q2237-I-0305 
(0stensen et al. 1995), and our current focus is on the po- 
tential for rapid intrinsic variability to cause a false trigger, 
we calculate model intrinsic light-curves with power spectra 
P„ (X 1/'" over a broad range of time-scales: between 1 and 
1000 days. We consider a 2-D parameter space of intrinsic 
light-curve variance (a|) and power-spectrum index 7/, and 
look for regions of this space whose values, when combined 



with microlensing models produce variability statistics that 
are inconsistent with observation. We first describe a method 
for computing the likely-hood (given assumed values for ai 
and 71) for a statistic / defined such that larger values of / 
imply a large intrinsic contribution to variability. The likely- 
hood is then compared to the observed value f{obs) for that 
statistic. We place limits on the values for erf and 7/ using 
4 different, but related light-curve statistics /. 

Pairs of single image light-curves can be combined to 
produce both difference and additive light-curves. While in- 
trinsic and microlensed variability cannot be directly dis- 
tinguished, intrinsic variation is only present in the latter. 
In addition, microlensed variability in different images is 
independent. Variability statistics for difference and addi- 
tive light-curves in the absence of intrinsic source variability 
should therefore be identical (on average). Thus comparison 
with the data of variability statistics for these two curves, 
averaged over 6 image pairs (only three are independent) 
provides a means of discrimination between different mod- 
els for intrinsic source variability. 

The statistic / is calculated for additive {/a) and dif- 
ference (/d) light-curves. Light-curves that include both in- 
trinsic and microlensed fluctuation exhibit higher variabil- 
ity in the additive light-curve (/a > /d). By comparing 
fA-D = fA ^ fo calculated from the monitoring data with 
values obtained from mock observations for different as- 
sumptions for the rate and amplitude of intrinsic variability 
(a'j, 7/), we compute the probability of obtaining the ob- 
served fA-D{obs) as a function of aj and 71. Microlensing 
statistics vary with effective transverse velocity and sam- 
pling rate. We have therefore computed this probability as 
a function of effective transverse velocity, using mock data 
sets computed from light-curves with the observed sampling 
rate: 

P{}a-d < fA-Diobs)\ai,ji,Veff,{m)) (1) 

We present calculations using mock observations that in- 
clude components of systematic {a^ya = 0.04) and ran- 
dom {ass) photometric uncertainty. To obtain P{fA-D < 
fA-D{obs)\ai,yi) Eqn. |l] is convolved with p„(t)e//|(m)): 

PUa^d < fA-D{obs)\ai,'yi) 

= J P{fA-D < fA~Diobs)\ai,'yi,v^ff, (m)) (2) 
xp„(veff\{m))dv,^ff 
Thus, the probability that 7/ < 70 and ai < ai,o is 
P('yi < 7o & CT/ < ai^olfA-oiobs)) 

= J P{fA-D < fA-D{obs)\ai,"/i)pprior{'yi,(^i)d"/i dai (3) 

We assume that Pprior(7/, o"/) oc 1. 

Fig. |l] shows plots of the 1, 5, 10, 90, 95 and 99% con- 
tours for Eqn. ^ calculated for 4 different variability statistics 
/, which we describe below. Statistics i), ii) and in) were 
calculated using the full set of observations. Statistic Hi) 
employs a derivative analysis and was therefore calculated 
from a modifled data set in which some adjacent points were 
averaged to reduce noise. We note that this reduces the po- 
tential of statistic Hi) to measure intrinsic variability of high 
frequency. 

We include recent data from the OGLE collaboration 
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a) Limits from 
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Figure 1. Exclusion plots for the light-curve variance and power-spectrum index 7/ for intrinsic variability in Q2237-t-0305. Limits are 
imposed by consideration of a) light-curve variance, b) light-curve variance about local mean, c) histogram of light-curve derivatives, 
and d) autocorrelation time-scales. 



(Wozniak et al. 2000a, b; OGLE web-page), which has a 
much denser sampling rate (~1 point per week). This in- 
creases the sensitivity to measuring high frequency intrinsic 
variability. While the OGLE data is in V-band, data taken 
prior to 1996 is in the r and R-bands. Microlensing is colour 
dependent when the source is near a caustic (Wambsganss 
& Paczynski 1991). In Q2237-I-0305 (5" < -qo), this is only 
the case a small fraction of the time and we therefore as- 
sume that the microlensing statistics are consistent between 
V and R-bands. Any colour dependence of intrinsic variabil- 
ity is expressed in the additive light-curve. However, Giveon 
et al. (1999) find that the rms amplitudes for the sample 
of Palomar-Green quasars for B and R bands differ by only 
~ 15%. We therefore assume that the power-spectrum pa- 
rameters we obtain apply equally well to intrinsic source 
variability in both the R and V-bands, and use all the avail- 
able data points without application of a colour correction. 

i) Light-curve variances: In the absence of intrinsic variabil- 
ity, the variance of additive (V^) and difference {V^) light- 
curves are identical. However, the variance of light-curves 



that include both intrinsic (erf) and microlensed fluctuation 
obey Va ~ + 2a j . Therefore the statistic 

f = VA~VD (4) 

can be used to place limits on intrinsic variability. From Fig. 
|l^ we find that at the 95% level the variance of intrinsic vari- 
ation is limited to be aj < 0.28 (for all 7/ considered). The 
light-curve variance does not discriminate between different 
values of 7/. 

n) Local light-curve variance; To probe the shorter fre- 
quency intrinsic variation, we have modified statistic i) in 
the following way. Each light-curve is defined by Nobs ob- 
servations, with ti and AMi representing the time of ob- 
servation and light-curve magnitude respectively. For each 
light-curve we calculate the variance about a local mean de- 
fined by: 

{AM),ocal(t) = ^ (5) 

where 
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n{u ~ t) 



1 for - 185 < ti - t < 185 days 
otherwise 



(6) 



The variance V"^(iomZ) — {{AMi — {AAl)iocai{t))^) measures 
fluctuations in the hght-curve with durations ^ 1 year. Fig. 
|l|3 is a plot of 90, 95 and 99% contours from Eqn. ^ calcu- 
lated for 



/ = VA{local) — Voilocal) 



(7) 



Tighter limits are imposed on the amplitudes of more rapid 
variations. As an example, the combination of 77 = —1 and 
ai — 0.2 is ruled out at the 95% level. 

in) Light-curve derivatives: Light-curve derivatives can be 
used to study quasar variability. Histograms of light-curve 
derivatives that result from difference and additive light- 
curves have a shape that differs more where the intrinsic 
variability is higher. R is defined as the value required to 
multiply the difference light-curve derivatives such that the 
KS difference between the resulting histogram and the addi- 
tive light-curve histogram is minimised. We have computed 
Eqn. m with 



f = R, 



(8) 



and plotted the 5, 10, 90, 95 and 99% contours in Fig 
For intrinsic variability having a variance that is small with 
respect to that of microlensing the statistic R demonstrates 
similar behaviour to (local). However, since the distri- 
bution of microlensed light-curve derivatives has a different 
shape to the distribution resulting from a power spectrum, 
the minimised KS difference is large if the intrinsic variance 
is high. This method is therefore inapplicable if the intrinsic 
variability is comparable to the microlensed variability. On 
the other hand, if intrinsic variability is assumed to have a 
small amplitude with respect to microlensed variability, R 
provides the tightest limits on crj. 

iv) Autocorrelation time-scales: Variability can also be stud- 
ied through computation of the autocorrelation function 
(ACF): 



ACF(r) = {(AMit) - {M)).{AM{t + r) - (A/»>. 



(9) 



We estimate the ACF by calculating the z-transformed dis- 
crete autocorrelation function (Alexander 1997). The width 
of the ACF peak (which we take to be the zero crossing time) 
can be used to define a correlation time-scale (T). Signifi- 
cant intrinsic variability will effect T. The 4th statistic is the 
difference between the correlation time-scale for the additive 
(Ta) and difference (Td) light-curves 

} = Ta-Td. (10) 

The 50% and 90% contours of Eqn. ^ are shown in Fig |^. 
The autocorrelation time-scale suggests that the amplitude 
of intrinsic variability is small. However, it does not place 
useful limits on ai or 7/ since due to its larger variance, 
the auto-correlation time-scale is dominated by microlensing 
variability. 

These results suggest that the intrinsic variability has 
small amplitudes and long durations. 
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Figure 2. Correlation plot of light-curve derivatives for a model 
light-curve of image A. The plot shows the cases of a full 2- 
dimensional convolution of the source intensity profile with the 
microlensing magnification grid vs. the single source-line approx- 
imation. The solid and dashed lines show the mean and iblo" 
correlations. The light line is the line of equivalence. 



4 DETERMINATION OF A TRIGGER FOR 
TOO OBSERVATIONS OF AN HME 

There are many properties of the Q2237-f 0305 light-curves 
that could be used as a trigger, ie. features in the observa- 
tions that signal the onset of an event and provide informa- 
tion on when that event occur. These range from a very sim- 
ple property such as the light-curve derivative, to a complex 
procedure involving fitting functional forms or even sections 
of model light-curve to the monitoring data. In the instance 
of perfect monitoring data, a simple trigger would provide 
ambiguous results, while more complex procedures will be 
more accurate. However, given the finite accuracy of mon- 
itoring data, simple procedures will provide more definite 
predictions since photometric noise is amplified in higher 
derivatives. We therefore investigate the light-curve deriva- 
tive as a tool to predict caustic crossings. 

There are two important quantities: i) What is the light 
curve derivative that signals the onset of a caustic crossing, 
and ii) given that derivative, what is the separation between 
the current observation and the event peak. We calculate 
probability functions for the separation between the obser- 
vation of a rapid light-curve rise and the subsequent HME. 
Note that while it would still be interesting, a drop in ob- 
served flux indicates that the HME has already peaked. 

We define a trigger as a 3-point light-curve derivative 
that lies above a value T . Due to the sparse sampling of 
monitoring campaigns, higher order numerical derivatives 
may have a baseline that is large with respect to the time 
between the observation and a subsequent event. To inves- 
tigate the information on future HMEs contained in trigger 
using model light-curves M{t), we step along the light-curve, 
looking for a derivative M > T, and examine the subsequent 
light-curve for an HME. Upon finding a trigger, our algo- 
rithm locates the next event and resumes the search from 
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there. Thus each event is only found once during a pass along 
the light-curve. 

In the calculations that follow, caustic crossing HMEs 
are defined to occur at the time when the caustic crosses 
the centre of the source. We note that this is not the same 
as the time of the maximum event amplitude, but it is the 
more relevant time with respect to spectroscopic observa- 
tions. However these points differ by only a few days. A cusp 
event is defined as occurring when the light-curve reaches its 
local maximum. 

We have assumed a Gaussian source profile. Because it 
is impractical to compute the full integrations of extended 
sources over the range of source size and microlens mass 
considered, we have approximated the extended source light 
curve by an integral along a single source line (eg. Witt 
& Mao 1994). To check the validity of this approximation 
we have computed light curves for a Gaussian source using 
both the 1-D approximation and a full 2-D integration over 
a 50x50 grid as described in Wyithe & Webster (1999). The 
light-curves were computed for a 2x 10^^ cm source, which is 
the upper limit of source size for Q2237-I-0305 (Wambsganss, 
Paczynski & Schneider 1990; WWTMOO), and microlenses 
of IM0 . Fig. ^ is a plot of light-curve derivatives calculated 
from the 1-D approximation vs. those calculated from the 2- 
D convolution (the mean and ±la lines are shown). The plot 
shows that there is some scatter, but that the approximation 
is a reasonable one at all gradients. The plot demonstrates 
the validity of the approximation in the case of the upper 
limit for source size (where it is least valid). The relationship 
is tighter fo r sm aller sources. 

(limiting our attention to images A and C) 



In Sec. 4.1 



we discuss the calculation of a function that describes the 
onset of a HME based on observed light-curve derivatives. 
Sees. <L2, 13 and discuss the systematic effects of un- 
knowns on the accuracy of the function obtained, and Sec. 
discusses t he a ppropria te sa mpling rate for monitoring 



4.5 



data. In Sees. 4.6, 4.7 and 4.8 



we present triggering func- 
tions specific to different classes of HME, discuss the opti- 
mal trigger and briefly discuss a method for determining the 
function appropriate for an isolated light-curve derivative. 



4.1 A general function for HME triggering 

We would like to know the trigger that will best anticipate 
an oncoming HME. We consider a grid of triggers T and pe- 
riods of time P following the trigger. For each trigger T we 
calculate the probability of observing a HME within each 
period P for a collection of assumed effective galactic trans- 
verse velocities, source sizes and mean microlens masses: 



Ft (PIT, S, {m),v,ff] 



(11) 



This function is calculated for HMEs located on 1000(m)r;o 
of light-curve, and is evaluated for combinations of 50 val- 
ues for Veff {Okmsec~^ < Veff < 2000 km sec~^, 20 val- 
ues for (m) (O.OOIM0 < (m) < IOOMq) and 20 values for 
S (10"cm < 5" < 10^® cm). The probability for observing 
a HME in the period P following a trigger T can then be 
obtained: 



Ft {P\T) 



dm dS dveff {ps iS\{m),Veff) Pm i{m)) 



o 



O 




5 10 15 

Trigger (Mags/year) 

Figure 4. Ft(P\T) (left) and dFT{P\T)/dP (right) for combi- 
nations of the different assumptions used in the calculations of 
p„, ps and pm. The solid and dashed lines correspond to small 
(SE) and large (LE) errors respectively. The light and dark lines 
correspond to 7^ > and 7^ < 0. The microlensing parameters 
correspond to those of image A. 



X p„ («e//Km)) Ft {P\T,S, {m),v,ff)) (12) 

The upper and lower plots of Fig. ^ show contours of 
Ft{P\T) (left) and dFT{P\T)/dP (right) for images A and 
C respectively. The sampling rate was 1 point per 2 days. 
The functions quantify the intuitive notion that steeper ob- 
served light-curve gradients imply a more imminent event. 
In addition, the plots show that the separation of the trigger 
and event peak is less sensitive to trigger size for larger trig- 
gers. There is little difference between the functions Ft{P\T) 
calculated for images A and C. 



4.2 The effect of model assumptions on the 
triggering function 

We have investigated the systematic dependence of Ft{P\T) 
on our different model assumptions. Fig. 4 shows Ft{P\T) 
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Figure 3. Ft{P\T) (left) and dFT{P\T) / dP (right). The microlensing parameters correspond to those of image A (top) and C (bottom). 



for combinations of shear direction, smooth matter percent- 
age and the photometric error in pubhshed monitoring data 
used to estimate Pm, Ps and Pv Plots are shown for im- 
age A, and were computed with a samphng rate of 1 point 
per 2 days. These plots show that the triggering function 
is quite insensitive to these systematic uncertainties. This 
is not surprising since the effective transverse velocity was 
determined for these different cases from a derivative analy- 
sis. In the remainder of this paper we consider only the case 
where 7a > 0, ac —0% and the photometric errors are small 
(SE). 

4.3 The effect of intrinsic source variation on the 
triggering function 

The gradient of a light-curve is influenced by intrinsic vari- 
ation of the source. We include intrinsic source variation 



in our simulations according to a power-law shaped power- 
spectrum with an index of 7 = —2 (following the findings of 
Giveon et al. (1999)) and periods between 1 and 1000 days. 
Sec. suggests that an intrinsic average light-curve variance 
of ai — 0.2 is higher than expected and that ctj = 0.4 is 
very unlikely. Fig. ^ shows Ft{P\T) calculated with intrin- 
sic fluctuation included for the cases of aj — 0.2 (top), and 
ai = 0.4 (bottom). Contours for Ft{P\T) with no intrinsic 
variation are included for comparison (dashed lines). The 
former has minimal effect, while a large intrinsic variance 
causes false triggers (no HME) at lower trigger values. We 
conclude that intrinsic variation should not have an impor- 
tant observable effect on triggering of HME observations in 
Q2237+0305. 
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Figure 6. Contour plots for Ft{P\T) computed from model light-curves having a size distributed according to our best estimate S 
(left), S" = 2S and S" = 55. There is no smooth matter in these models, small photometric errors (SE) and 7Ai7fl > are assumed. 



4.4 The effect of an incorrect source size on the 
triggering function 

The determination of probability for the quantity Vef / (m) 
from the Q2237-I-0305 monitoring data is quite robust. How- 
ever, the source size estimate is more uncertain, having been 
derived from a single poorly sampled HME. The small num- 
ber of observations describing the 1988 peak (Irwin et al. 
1989; Corrigan et al. 1991) introduces the potential for a 
systematic error in the estimate of source size equal to the 
ratio of true event length divided by the inferred event length 
of 52 days (twice the estimated rise time). This can be com- 
pared to the ~100 day separation of the two points that 
provide an upper-bound to the event duration. This system- 
atic error in the estimate of source size is therefore smaller 
than a factor of 2. Fig. ^ shows the effect (on F) of replacing 
Ps {S\{m),Veff) with 

p's{S'\{m),v,ff)dS' = ^Ps{2S\{m),v,ff)dS (13) 
(centre) and 

pi {S'\{m),v,ff)dS' = ^ps{5S\{m),v,ff)dS (14) 

(right). Contours of the resulting functions Ft{P\T) are 
shown together with Ft{P\T) computed from our best esti- 
mate of Ps for comparison (left hand figure). Inspection of 
Fig. 1^ shows that if the systematic error in the estimation of 
source size is the worst case of a factor of 2, then the fore- 
casting of events will be effected if the trigger is greater than 
~ 5 magnitudes per year. However, the forecasting of events 
is far more uncertain for all triggers if the estimate of source 
size is incorrect by a factor of 5. While the larger source 
produces flatter peaks, longer events and hence a smaller 
maximum derivative, the large scale derivatives occur earlier 
due to parts of the source having moved into close proxim- 
ity of the caustic sooner. Hence the larger assumption for 
source size predicts a larger separation between trigger and 
corresponding event peak. 



4.5 What is the appropriate sampling rate for 
monitoring of Q2237+0305? 

It would be ideal to monitor the Q2237-I-0305 daily. However 
as this is not always possible we investigate the success of 
predicting the onset of a caustic crossing using different sim- 
ulated sampling rates. The light-curve derivatives have been 
calculated with spacings of 2, 8, 16 and 32 days to mimic a 
trigger obtained from regular monitoring. Fig. ^ shows plots 
of the functions Ft{P\T) computed with the above sampling 
rates. Here we have defined the delay of the event following 
the trigger with respect to the central (of the three) obser- 
vation used for the calculation of the derivative. We find 
that trigger probabilities calculated with spacings of up to 
~ 16 days provide similar predictions. However as the spac- 
ing is increased, predictions diverge at the large triggers. In 
particular, for the larger sampling rates, very large triggers 
are rarely observed because the sampling length is larger 
than or comparable to the event time-scale. The functions 
Ft{P\T) are only useful for P larger than the sampling rate. 
Note that in realistic situations the sampling rate will not 
be regular due to a variety of problems (weather etc.). How- 
ever Ft{P\T) predicts light-curve behaviour based upon the 
observation of a single derivative. A regular sampling rate 
over an extended period is therefore not essential. 



4.6 Triggering functions for different classes of 
HMEs 

Due to the appearance {+ve event) or disappearance {—ve 
event) of a pair of critical images, caustic crossing HMEs are 
asymmetric in time. In particular, the light-curve remains 
at a higher level following the event if a new pair of critical 
microimages is created. In addition, the creation of critical 
images causes the light-curve to rise much more steeply than 
it falls. Where there is some knowledge about the type of an 
impending event (ie. whether the critical images are appear- 
ing or disappearing) the triggering function can be deter- 
mined more accurately. A third possibility is that an HME 
may not involve the source crossing a caustic but be due to 
the source moving past a cusp. This type of event is not in- 
teresting for a spectroscopic measurements (eg. Wambsganss 
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Figure 8. Contour plots for Ft{P\T) (left) and dFT(P\T)/dP (right) for images A and C. Ft(P|T) and dFT(P\T)/dP are shown for 
3 separate classes of events: cusp (top), caustic crossings with disappearing (centre) and appearing (bottom) critical images respectively. 
There is no smooth matter in these models, small photometric errors (SE) and 7A)7fl > are assumed. 



& Paczynski 1991), and so forewarning would be useful. We 
have determined the functions Ft^_ (P|r), dFT+{P\T) /dP; 
Ft_ {P\T), dFT_ {P\T)/dP; and Ftc{P\T), dFTc{P\T)/dP 
as before, looking for triggers that portend an event due to 
only appearing {+ve) or disappearing {—ve) critical images, 
or a cusp respectively within a given period. Note that 

Ft+{P\T) + Ft_ {P\T) + Ftc{P\T) = Ft{P\T), (15) 

so at large P, Ft_^ , Ft^ and Ftc describe the relative frac- 
tions of different event types found by that trigger. 

Fig. I shows contours of fT^(P|r), dFT+{P\T) /dP; 
Ft_ (P|r), dFT_ {P\T)/dP; and Ftc {P\T), OFtc {P\T)/dP. 
Because the rising side of a — «e event (or a cusp event) is 
not as steep as that of a +ve event, very large triggers are 
not likely to precede a —ve event (or a cusp event). How- 
ever, in the case where a large trigger is observed prior to 
a —ve event, the delay P is small. In contrast, a small trig- 
ger portends +ve and —ve events in approximately equal 
numbers. Similarly dFT_{P\T)/dP has a maximum for an 
event arriving ~ 3 weeks after a trigger of ~ 6 magnitudes 
per year, however dFT^{P\T)/dT is larger at higher T and 
smaller P. In addition, following a trigger ^ 5 magnitudes 
per year, +ve caustic crossing events are more likely, and 
are expected sooner than —ve caustic crossing events. 

FTciP\T) and dFTc{P\T)/dP behave similarly to 
the corresponding functions for —ve caustic crossings. For 
smaller triggers, the delay until an event is larger. The simi- 



larity in shape of Ftc {P\T) and Ft_ (P|T) will make distin- 
guishing between a cusp and —ve caustic crossing difficult 
prior to the event peak. On the other hand Ftc{P\T) < 
Ft_ {P\T) and so cusp events are less common. 



4.7 Selecting a trigger 

While a derivative observed from monitoring data can be 
analysed as a potential trigger in real time it is useful to 
have a pre-conception of the optimal trigger. In particular, 
a program could demand that the trigger be larger than the 
maximum derivative expected for cusps. However, if this ap- 
proach were taken, then some —ve caustic events may peak 
before a TOO was triggered. A large trigger will always yield 
the desired caustic crossing, however HMEs may be missed 
before such a trigger is observed (and HMEs are rare). On 
the other hand if the TOO is triggered on every low deriva- 
tive then the success rate will be lower, but the observation 
of the next HME is guaranteed. A balance must be struck 
between the quantity of available TOO observing time, and 
how long one is willing to wait to obtain the observations. 

We have calculated the cumulative probability of ob- 
serving a maximum light-curve derivative prior to the peak 
maximum: 

Pj^,{M<Mo) = 
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Figure 5. Contour plots for Ft{P\T) computed from model light- 
curves which include an intrinsic component of variability with 
a power-spectrum having a power-law form with index 77 = —2 
and a variance crj = 0.2 (top) and aj = 0.4 (bottom). The dashed 
contours in each case show _Fy(P|T) in the absence of intrinsic 
variability. There is no smooth matter in these models. The mi- 
crolensing parameters correspond to those of image A (7^1 > 0), 
and the photometric errors were assumed to be small (SE) . 
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Figure 7. Contour plots for Ft{P\T) computed from model light- 
curves with sampling rates of 2, 8, 16 and 32 days. There is no 
smooth matter in these models, small photometric errors (SE) 
and 7^1, 7^ > are assumed. 
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Figure 9. Plots of the probability for the maximum derivative 
preceding each a light-curve peak. Functions are given corre- 
sponding to derivatives that precede cusps (solid line) , and caustic 
crossings with disappearing (dotted line) and appearing (dashed 
line) critical images. Only peaks larger than 0.1 magnitudes above 
the local minima on both sides are considered. 



jdmjdsjd..,, (p.(5|M,..,,)p„(M) 

X p„ («e//!(m))P^y (M < Mo\S, {m),«e//)) (16) 

using a sampling rate of 1-week. We have only presented the 
case of models with no smooth matter, small photometric er- 
rors (SE), 7a, 7s > and no intrinsic source variability. We 
have looked at the probability P for cusps, Pj^j^{M < Mo); 
—ve caustic crossings, Pj^j {M < Mo), and +ve caustic 
crossings, Pmj^{M < Mo)- Fig. |^ shows these functions; 
solid, dotted and dashed lines correspond to Pm^ {M < Mo), 
Pm_{M < Mo), and Pm^{M < Mo) respectively. 

Fig. ^ shows that a trigger of ~ 5 magnitudes per year 
will successfully trigger observations for most —ve caustic 



crossings, all +ve caustic crossings, and miss ~ 90% of cusp 
events. The source size has been estimated from the 1988 
event in R-band. The trigger should therefore be from obser- 
vations made in R-band since then the triggering function 
is more accurately determined. Also, accretion disc mod- 
els predict that the source should be smaller in higher fre- 
quency bands: monitoring and triggering at a higher fre- 
quencies therefore provides more rapid variations, but gives 
less warning of the impending event. 

4.8 An alternative approach 

Recall that the method to determine F steps along the light- 
curve until a derivative is located that is greater than T, 
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the search is then resumed after the subsequent event. The 
probabiUties therefore refer to the fractions of events of each 
type that are located on a Ught-curve and have a derivative 
greater than T on their leading side (integrated over effective 
transverse velocity) . This calculation is most useful when an 
immediate history of image brightnesses is available, so that 
it is known that the current derivative is the first above T. 

Alternatively, the a-priori probability of observing an 
event following a light-curve trigger Ti AT can be computed 
by finding the derivative and separation of the derivative 
and corresponding event at all (equally spaced) points on 
the light-curve. This calculation does not take into account 
the recent history of the light-curve, and should be applied 
to an isolated observation (eg the first point of the observing 
season. The resulting Ft^{P\T±AT), Ft_{P\T±AT) and 
Ftc (-P|r± at) describe the likely-hood of the class of event 
following the trigger in the absence of previous light-curve 
information rather than the relative numbers of events of 
each type that are located on a light-curve. 



5 A MOCK TARGET OF OPPORTUNITY 
EXPERIMENT 

In this section we combine high resolution microlensing mag- 
nification grids with a thermal accretion disc source to pro- 
duce model microlensed light-curves in 11 bands. We then 
produce example target of opportunity spectroscopic obser- 
vations for hypothetical light-curve triggers. By looking at 
the microlensing induced change in spectral slope we esti- 
mate the chance of success for a three observation TOO pro- 
gram. The three spectra are taken at intervals fixed at the 
time of triggering such that they are likely to observe the 
light-curve peak (ie. according to F{P\T)). We emphasise 
that this section does not predict the colour changes that 
will be measured, rather it demonstrates that a physically 
motivated source with a spectrum that is radius dependent 
will produce colour change that could be detected by target 
of opportunity observations. 




Figure 10. Cross-sections of tiie intensity profiles in the U (ligiit 
line), R (thick dark line) and K (thin dark line) for the thin 
accretion disc source. 7^ is the intensity per unit frequency on the 
surface of the accretion-disc, normalised to the maximum value 
in U-band. The profile is shown for a disc viewed face on, hence 
the deficit of flux at the centre. 



Table 1. The frequencies at the centre of the bands for which 
source intensity profiles and light-curves have been calculated. 



u {Hz) 


band 


5.50x1015 


F160AW 


3.90x1015 


F218AW 


2.77x1015 


U (F300W) 


2.06x1015 


U'(F380W) 


1.73x1015 


B 


1.48x1015 


V 


1.20x1015 


R 


1.01x1015 


I 


6.47x101'* 


J 


4.90x10" 


H 


3.63x1014 


K 



5.1 The model accretion disc 

We use a standard thin accretion disc with constant accre- 
tion rate and no advection of heat. The energy generation 
per unit area as a function of radius is given by 

„ 3 GMbhM „ , , 

Q = - 3 Rnir). 17 

Mbh is the black hole mass, M is the accretion rate, r the 
radius and Rr is a correction factor. Rr combines outward 
advection of energy associated with the angular momentum 
flux and relativistic effects (Page & Thorne 1974, Krolik 
1998). To calculate the flux generation at the surface of an 
accretion disc, the relativistic effects of Doppler beaming of 
radiation, gravitational redshifts, and bending of photon tra- 
jectories must be considered. These are included through the 
calculation of a relativistic transfer function (Cunningham 
1975). The following assumptions are made to compute the 
transfer function: i) the accretion disc is thin (disc height 
much smaller than disc radius), ii) the gas follows prograde 
circular orbits outside the marginally stable radius rma , Hi) 
the disc is flat and lies in the equatorial plane of the black 



hole and iv) the gas emits isotropically in its rest frame. 
Note that this model does not emit inside Tma- The disc 
model is discussed in more detail in Agol & Krolik (1999). 

5.2 The model accretion disc light curves 

Several attempts have been made to understand the accre- 
tion disc in Q2237-)-0305 by comparing model light curves 
with observed microlensing events (eg. Ranch & Blandford 
1991; Jaroszynski et al. 1992; Jaroszynski & Marck 1994; 
Czerny et al. 1994). The approach of these papers is to dis- 
cuss the luminosity of the quasar considering the predicted 
magnification of the macro-images as well as the limits on 
source size imposed by microlensing observations. In this 
paper we compute light-curves for the accretion disc in 11 
frequency bands (see table The quoted frequency is the 
band centre in the quasar rest frame, the first 4 filters are 
HST filters. Fig. |l^ shows the cross-sections of the inten- 
sity profiles in the U (light line), R (thick dark line) and 
K (thin dark line). The source size (5) must be interpreted 
in terms of its intensity profile. We define the source size 
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such that the event length produced by the source is equiv- 
alent to that produced by a top-hat profile of diameter 5*. 
We have used the example of a black hole mass of 10* M© 
which at 1.2x10^^ Hz (R-band) corresponds to a source ra- 
dius of 5 X 10^* cm. This model suffers the difflculties that 
were pointed out in the case of a thermal accretion disc by 
Ranch & Blandford (1991), namely that it severely under- 
estimates the observed flux. However our focus here is to 
combine a detailed accretion disc model with high resolu- 
tion microlensing models to simulate microlensing induced 
spectral change during a caustic crossing. 

The source profile is computed over a nested grid of 
cells (with side length SOOrg) down to a resolution of 0.5 
Tg in each of the 11 bands. The grid contains 4864 cells 
which lie on 192 parallel source tracks. The light-curves were 
computed using the method described by Wyithe & Webster 
(1999), and computed at a resolution of 1.5xl0~''7;o along 
the source line. At an effective galactic transverse velocity 
of 400 fcm see" ^, this corresponds to a temporal resolution 
of ^ 0.5 days. 



5.3 Examples of model HME observations 

We calculate the microlensed light-curves of an accretion 
disc source having a central black-hole mass of 10* M© 
{Sr ~ 10^^ cm in R-band). The mean microlens mass in 
the simulation was (m) = O.IMq and the microlensing pa- 
rameters corresponded to image A (7^ = 0.4). The results 
for S and (m) of WWTMOO suggest that these parameters 
are respectively near the upper and lower limits expected for 
Q2237+0305. The colour changes found are therefore con- 
servative for this class of disc model. Example light-curves 
with the aforementioned parameters are presented in Fig. 
Three curves are shown corresponding to the source profiles 
in Fig. |l^. The light-curves show magnitude change plus an 
offset equal to the difference in unlensed flux. 

Ensembles of mock experiments are produced for two 
hypothetical triggers: T = 5 and T = 15 magnitudes per 
year. Below we simulate a TOO comprising 3 spectroscopic 
observations taken at intervals designated at the time of the 
trigger. The first two observations are scheduled according 
to when the triggering function dFriT, P)/dP predicts an 
HME: 2 — 4 weeks following the small trigger and 1 — 2 weeks 
following the large trigger. A third observation is made 10 
weeks following the trigger to provide a comparison with 
the spectrum of a (hopefully) non-differentially magnified 
source. 

We assume a likely effective galactic transverse velocity 
of 400 fcm see" ^, and trigger observations from the R-band 
light curve in Fig. |l^. Examples of the spectrophotometric 
simulations are presented in Figs. ^ and In each case, 
four spectra are shown in the upper panel, three microlensed 
spectra as well as the unlensed spectrum (light-line) for com- 
parison. The lower panels show close-ups of the events in the 
U, R and K bands, the vertical lines show the position of 
the three observations on the light-curve and have line-styles 
that correspond with the spectra in the upper panel. The po- 
sition of the trigger is denoted on the R-band light-curve by 
a dot. 

In the example shown the black hole mass of lO^Mo 
produces a larger source than microlensing observations pre- 
dict is most likely. The triggering function, which is calcu- 
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Figure 13. Examples spectrophotometric simulations that result 
from a trigger of T = 15 magnitudes per year from the R-band 
light-curve shown in Fig. Presentation as per Fig. ^ 
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Figure 11. Model light-curves of image A for the U-band (light line), R-band (thick dark line) and K-band (thin dark line) sources. 
The microlensing model had a mean microlens mass of (m) = O.IMq, and the source black-hole mass was Mgjj = 10^ Mq. 



lated using the most likely source size therefore underes- 
timates the delay between the trigger and event peak in 
several cases (particularly for the larger trigger) because 
larger sources experience HME scale derivatives earlier than 
small sources. However, the spectral change is at an ob- 
servable level in nearly all cases. The size of the spectral 
change observed during the HME is qualitatively dependent 
on whether the event is a +ve or a —ve caustic crossing. 
In the —ve case, the first two observations are taken during 
the event. The final observation is made following the dis- 
appearance of the critical images at a time when the source 
has little differential magnification. During the first two ob- 
servations, the higher frequencies (emitted from smaller re- 
gions) are preferentially magnified, and so the spectrum 
is stronger at the higher frequencies. In contrast, the +ve 
events have initial observations that are taken before the 
event has peaked. There are no critical images of the inner 
regions and hence they are not as strongly magnified. The fi- 
nal observation has critical images of the entire source, and 
so an equal magnification in all bands. In both cases, the 
spectra plus a well sampled light-curve in a single band will 
allow determination of the relative sizes of the emitting re- 
gions at all frequencies which are in contact with the caustic 
for the first two observations. The trigger of T = 5 magni- 
tudes per year finds +ve and —ve events in approximately 
equal numbers. However, the trigger of T = 15 magnitudes 
per year finds only +ve events because —ve events do not 
reach such large light-curve gradients before the event peak. 

To quantify the effect of the model TOO experiment, 
we calculate the distribution of spectral slopes for each of 
the three observations in an ensemble of model TOOs. The 
triggering function was applied to 50r?o (~ 500 years at 
400 km sec~^) of light- curves for the source and microlensing 
model discussed above, and continuum spectra calculated for 
three observations following the trigger. The 5 magnitudes 
per year trigger precedes 159 HMEs in these light-curves. 
However the larger trigger (15 magnitudes per year) only 
precedes 80 HMEs since such a large light-curve derivative 
is unlikely to precede —ve caustic crossings or cusp events. 



We have calculated the change in spectral slope a (over the 
wavelength range F160AW — U') for each event, which is 
defined as: 



log (F^FimAW)) - log [F^jU')) 
log {v{FimAW)) - log {u{U')) 



(18) 



Fig. |l^ shows the distributions of spectral slope for the three 
observations following triggers of T = 5 (upper panels) and 
r = 15 (lower panels) magnitudes per year. The results ob- 
tained for the two triggers are reassuringly similar. There is a 
large spread about the unlensed a = 1 in the spectral indices 
obtained from the first two observations. As mentioned, dif- 
ferential magnification from microlensing can produce both 
an increase and a decrease in the spectral slope. This is 
refiected in the distributions which are approximately sym- 
metric about a = 1. The third observation at 10 weeks is 
intended as a control and Fig. |l4| shows that the source 
is rarely subject to strong differential magnification at this 
time. The range of spectral indices for this observation is 
small, and its mean is consistent with the unlensed value of 
Q 1. 

Fig. ^ shows histograms of the larger difference (in ab- 
solute value) between observation 1 /2 and observation 3 fol- 
lowing triggers T = 5 and T = 15 magnitudes per year. For 
a successful measurement of microlensing induced change in 
spectral slope to be made, a change in a will need to be ob- 
served both at a level above observational uncertainty, and 
the expected spread of values for observation 3. Compari- 
son of Figs. |l| and |l| shows that this is true in most cases. 
To quantify the chance of success for the experiment we de- 
fine a TOO program to be successful when one of the initial 
two observations lies outside the variance of the third obser- 
vation (computed from the entire ensemble). We find that 
this definition yields encouraging success rates of ~ 80% and 
~ 88% for triggers of T = 5 and T = 15 magnitudes per year 
respectively. Of course if more than three observations were 
available, the chances of success would be enhanced. Unsuc- 
cessful TOOs may be due to triggering on a cusp event (more 
likely for T — 5, hence the lower success rate) or to an un- 
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Figure 12. Examples spectrophotomotric simulations that result from a trigger of T = 5 magnitudes per year from the R-band light- 
curve shown in Fig. In each case, the three spectra (flux per unit frequency vs. frequency) are shown in the upper panel. An unlensed 
spectrum is shown for comparison (light line). The lower panels show close-ups of the event in the U, R and K bands, the vertical lines 
show the position of the three observations on the light-curve and have line-styles that correspond with the spectra in the upper panel. 
The position of the trigger is denoted on the R-band light-curve by a dot. 
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Figure 14. Histograms of the microlensed spectral slope between HST bands F380W (2.06X lO^'^Hz) and F160AW (5.50x lO^^Hz). 
Histograms for three observations are shown. Top: 2, 4 and 10 weeks following an R-band trigger of 5 magnitudes per year. Bottom: 1, 
2 and 10 weeks following an R-band trigger of 15 magnitudes per year. 



usual caustic crossing whose peak lay outside the expected 
range. 

The above calculation considers only the effect of mi- 
crolensing on the spectral slope at high-frequency. While 
this is the regime where microlensing is expected to have 
the largest effect it is also interesting to see if similar infor- 
mation can be obtained at longer wavelengths, specifically 
from the spectral slope between the K and B bands: 

log{F4K))~log{F4B)) 
- log{u{K))~log{u{B)) ' ^ 

In this case the source is much larger and so the control ob- 
servation must be taken at a much later time. Fig. |l^ shows 
the histogram of spectral indices for three observations at 
2, 10 and 30 weeks following a trigger of T = 5 magnitudes 
per year. The control observation has a large spread (3 times 



that for the high frequency index), indicating that even after 
30 weeks, there is still some differential microlensing due to 
the caustic. The TOO experiment is much more likely to be 
successful at higher frequencies. 



6 CONCLUSION 

We have calculated the probability (triggering function) of 
observing a caustic crossing within a given period following 
the observation of a rise (trigger) in an image light-curve 
of Q2237-I-0305. We find that there is a well defined mode 
for this distribution, that the mode occurs earlier, and that 
the distribution peak becomes narrower for observations of 
larger light-curve derivatives. Thus, the more rapid the ob- 
served rise in monitoring data, the sooner a caustic crossing 
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Figure 15. Histograms of the larger difference between the microlensed spectral slope of observation 1 or 2 and observation 3. the 
spectral slope was calculated between between HST bands F380W (2.06x1015Hz) and F160AW (5.50x IQi^Hz) following an R-band 
trigger of Left: 5 magnitudes per year, and Right: 15 magnitudes per year. 
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Figure 16. Histograms of the microlensed spectral slope between K and B bands (3.63x lO^^Hz) and (1.73x lO^^Hz) for the three 
observations at 2, 10 and 30 weeks following an R-band trigger of 5 magnitudes per year. 



is expected, and the more accurate its predicted arrival time 
will be. Predictions of caustic crossings are aided by knowl- 
edge of whether a pending event will be due to a pair of 
appearing or disappearing critical images. In particular, an 
event that is due to appearing critical images is expected 
sooner than one due to disappearing critical images. The 
probabilities for cusp events are similar to those for events 
due to disappearing critical images, making distinguishing 



these difficult before the event event has reached its maxi- 
mum. However cusp events are far rarer than caustic cross- 
ings. 

We have checked the effect on our triggering function of 
different combinations of assumptions for the size of photo- 
metric error in the monitoring data, the fraction of smooth 
matter content and the direction of the galactic transverse 
velocity. We find little systematic dependence on these quan- 
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titles. We have placed limits on the Intrinsic variation In 
Q2237+0305 and shown that (except for the smaller trig- 
gers) the predictions should not be effected by Intrinsic 
source variation. A systematic underestimation (the largest 
possible) of x2 In our previous estimation of source size 
would produce an underestimate In the arrival time. The 
underestimation Is most serious for large triggers preceding 
a caustic crossing high magnification event. Predictions from 
smaller derivatives are less affected. 

Computation of the distributions of the separation be- 
tween a trigger and subsequent high magnification event 
for different sampling rates determines the appropriate sam- 
pling rate for the prediction of caustic crossings from regular 
monitoring. We find that the triggering function is insensi- 
tive to sampling rate for sampling rates of ~ 2 weeks or less. 
If caustic crossing probabilities are calculated from monitor- 
ing data with a higher sampling rate, the derivative may be 
calculated over a longer baseline to reduce the photometric 
uncertainty. 

We have computed the light-curves of a thermal accre- 
tion disc about a central black-hole of 10*A/q in 11 bands. 
We look for hypothetical triggers on the model R-band light- 
curves of T = 5 and 15 magnitudes per year and compute 
two model spectra at times most likely to surround the peak 
(determined from the triggering function). A third spectrum 
is computed 10 weeks following the trigger and is Intended 
as a standard (not differentially magnified) . We find that the 
third observation is subject to very little differential magni- 
fication and the computed spectrum has a slope a (between 
2.06 X 10^^ Hz and 5.50 x lO^'^ Hz) that is barely changed 
from that of the unlensed source (a ~ 1). However the ini- 
tial two spectra have spectral slopes that differ from that of 
the unlensed source by between -0.5 and 1. 

We have calculated the chance of success of a target 
of opportunity program that contains 3 observations. The 
success rate is defined as the fraction of an ensemble of 
mock target of opportunity observations that have one of 
the initial two observations with a spectral slope outside 
the variance in the third (calculated over the entire ensem- 
ble). We find that the success rate should be high, ^ 80%, 
which demonstrates that a predefined triggering and TOO 
observational strategy can achieve the goal of observing 
spectrophotometrlc change during a mlcrolensing event in 
Q2237 l | 0305. 
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